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Abstract 

The Hamiltonian description of the self-consistent interaction between an electromagnetic plane- 
wave and a co-propagating beam of charged particles is considered. We show how the motion can be 
reduced to a one-dimensional Hamiltonian model (in a canonical setting) from the Vlasov-Maxwell 
Poisson brackets. The reduction to this paradigmatic Hamiltonian model is performed using a Lie 
algebraic formalism which allows us to remain Hamiltonian at each step of the derivation. 
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I. INTRODUCTION 



The interaction between electromagnetic fields and a beam of charged particles exhibit 
a dynamics which is nowadays not fully understood, mainly due to the complexity inherent 
to the infinite dimensional phase space. Although some attention has been devoted to 
controlling those systems, a preliminary step is to shed light on the dynamics by analysing 
for instance phase space structures and transport properties. In this ambitious framework, 
reduced models have proved to be very valuable for this task. The reduction is of course 
guided by the physics of a particular setting. For instance, a beam of ultra-relativistic 
electrons interacts with plane waves (whether these are external ones or produced by the 
beam moving into an array of magnets, called undulator) has shown to be a way for producing 
a coherent light source. In what follows, we consider such a beam of electrons moving inside 
an undulator which produces a static (but non-uniform) magnetic field. The acceleration 
produced by the external magnetic field makes them emit a synchrotron radiation, which 
is self-consistently interacting with the particles. Under some resonance condition, the 
intensity of this electromagnetic wave grows exponentially and then saturates. In order to 
capture this effect, a one-dimensional Hamiltonian model has been proposed. The reduced 
Hamiltonian describes the evolution of the position 9j and relative momentum pj (around a 
resonant value) of the j-th particle. The particles interact through a wave described by 
its intensity / and phase 0. It reads : 



where {Oj,Pj) and (</>,/) are canonical pairs of conjugate variables. 

In the Free Electron Laser (FEL) configuration, this model has been derived from 
the equations of motion of charged particles and Maxwell's equations. Some approximations 
were involved during the course of the derivation, and were guided by the physics of the 
device, and among them, a specific form for the radiated fields (as a plane wave), a one- 
dimensional reduction (obtained by reducing the dynamics in the transverse plane), an 
expansion around a resonance fixed by the characteristics of the undulator. In addition, 
this Hamiltonian model was also proposed to describe the wave-particle interaction in other 
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In this article, we propose a derivation of the one- dimensional Hamiltonian ([T]) for the self- 
consistent interaction from the Vlasov- Maxwell equations in a Hamiltonian setting. This 
allows us to show that at each step of the derivation, the Hamiltonian structure of the 
problem is conserved. Another advantage is that the conserved quantity (namely the to- 
tal momentum of the system) is also easily deduced from the conserved quantity of the 
Vlasov-Maxwell equations. Our approach follows from this algebraic framework: Instead of 
working with the equations of motion for the derivation, we consider the Hamiltonian and 
its associated Poisson bracket (for an introduction, see Refs. . By using a canonical 

version, the Poisson bracket remains canonical (or generalized canonical in a broader sense). 
Therefore the main approximations and computations have to be done on a scalar function, 
the Hamiltonian which aims at simplifying the derivation. 

In Sec. HI] we recall some basics of the Hamiltonian formulation of the Vlasov-Maxwell 
equations for a continuous description of the particle distribution. In Sec. Illlt we apply 
to this Hamiltonian system the approximations and reduction necessary for the derivation 
of the reduced model ([T]) by expressing first the Vlasov-Maxwell system into a canonical 
setting (Sec. IIII Ap . then performing the one-dimensional reduction (Sec. IIIIBI) . changing 
the reference frame (Sec lIIICj) . and expanding the resulting Hamiltonian around a resonance 
condition (Sec. IIII DP . Finally, we deduce (Sec. IIIIEp a conserved quantity of the reduced 
model from a conserved quantity of the Vlasov-Maxwell equations by following the same 
procedure as for the reduction of the Hamiltonian. 



II. HAMILTONIAN FORMULATION OF VLASOV-MAXWELL EQUATIONS 

The dynamics of Hamiltonian ([T]) follows from Hamilton's equations for each pair of 
canonically conjugate variables. More generally the dynamics of an observable F (function 
of phase space coordinates {{9i,Pi},I,(f))) is given by : 

f-wn (2) 

where the Poisson bracket between two observables is given by 

In a continuous setting, this Hamiltonian model can be extended in a straightforward way. 
The beam is now described by a distribution function f{6,p) which constitutes a dynamical 



field, i.e., for eacli location in phase space {0,p), the density of particles f{0,p), labelled 
by the phase space coordinates of the particles, evolves dynamically. The one-dimensional 
Hamiltonian model is generalized from Hamiltonian ([T]) 
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H[f,i,<p] = 11 dedpfie, 



p) 



y + 2\/l sm 



(3) 



where the dynamical variables are now / and and a field of variables /(^, p). The dynamics 
of /, I and 9 are obtained using the canonical Poisson bracket: 



{F,G} = jj dedpfie,p) 

i.e. it leads to a Vlasov equation for / 



d_6£d_6G d_6£d_6G_ 
d^JfdeJf ~ deJfd^Jf 



+ 



dF dG OF dG 



dl 



dl 



(4) 



dt 



+ p^ - 2v7 coslt/ - 



1^ 
dp 



0, 



where we notice that df /dt = df /dt since we use an Eulerian description for the observables. 
This equation has been used to estimate quantitatively some features of the dynamics, like 
the derivation of a reduced dimensional model 8|] or the characteristics of the bunching in 
the saturated regime 

If / is a Klimontovitch distribution, that is the distribution function is a sum of Dirac 
representing some point particles: 



f(9,p) = j2sio-9,itmp-p,{t)), 



we recover the equations for Hamiltonian ([T]). 

The continuous formalism ©-(jl]) is particularly well-suited for an algebraic treatment 



of the dynamics (see e.g. III]). In what follows we use Vlasov-Maxwell equations to derive 
Hamiltonian ([3]). In order to do this we use a Hamiltonian formulation of these equations. 
First let us recall that the interaction between electromagnetic fields and charged particles 
(of normalized mass m = 1 and charge e = 1) is given as the sum of the kinetic energy of 



the particles plus the energy of the field [12 

H= 1 1 (iVp/(q,p)/T 



mm- 



d^q 



|E(q)P + |B(q) 



(5) 



where /(q, p) describes the distribution of particles in phase-space. Even though the kinetic 
energy of the particles and the electromagnetic energy appear to be decoupled in the Hamil- 
tonian, the interaction between the matter and the fields comes from the bracket which gives 



the dynamics : 



{F,G} 



+ 



+ / d^'q 



d^qd^p f 
d^qd^pfB ■ 
d^qd^p 



d_5F_ d_6G_d_5F_ d_6G_ 
dp 6f dq 6f dq 6f dp 6f 



d_5F_ d_6G 
dp 6f dp 6f 



S£dl SG_5Gdl 5F_ 
If dp ' Je ~ If dp ' Je 

6F\ 6G 6F 



Vx — ). — - — . (Vx — 
6B 6E 6E \ 6B 



(6) 



This bracket satisfies the antisymmetry property, the Leibnitz product rule and the Jacobi 
identity. Here, the Lie algebra on which this bracket operates is the set of smooth functional 
F[/(q, p), E(q), B(q)]. Using Hamiltonian ([5]) and the brackets IQ, Eq. ([2]) allows one to 
retrieve Maxwell's equations for E and B, as well as Vlasov equation for / : 



dt 



^f = {HJ} 



-V • V/ - (E + V X B) 
V X B - / c/V/, 



dp' 



dB 

'dt 



B = {H,B} = - V X E, 



where v is the velocity 



(7) 



We notice that the first line of the right hand side of Eq. refers only to the particles 
(and it is canonical), the second and third lines are the field-particle interaction terms (non- 
canonical terms) and the last line is a field-only term (which is also canonical). 

III. INTERACTION BETWEEN A PLANE WAVE AND A CO-PROPAGATING 
BEAM OF PARTICLES 



A. Expression of the Hamiltonian system in a canonical way 



The Vlasov-Maxwell equations (I5l)-(l6]) can also be described using the potentials instead 
of the fields 



13]. The Lie algebra is now a set of functional -F[/mom(q, p), A(q), Y(q)]. The 
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Hamiltonian and the bracket become : 



r r , r ivp i ly y a|2 

H = jj rfWmomVl + (p-A)2 + j d^q ^ ' ^'^ (8) 



d 6F d 6G d 6F d 6G 



dp 5 f mom dciSfmoni <9q 5/mom dp S fn 



f f6F 6G 6F 6G\ 
This can be obtained from Eqs. (jS])-® using the change of coordinates 

/(q,P) = /mom(q,P + A), 

E = -Y, 
B = V X A. 

We notice that this time, the Poisson bracket is canonical and there is no term in this 
bracket which couples the particles and the field. However, the coupling term is present in 
Hamiltonian ([8]). 

We translate the potential vector by a static A^(q), which is imposed externally (as 
in an undulator). We notice that a translation of A by a quantity A^ is a canonical 
transformation, which implies that the bracket is not changed. The new Hamiltonian 
reads : 

|Y|2 + 2V X ■ V X A + |V X A|2 

2 ' 

where we have dropped the constant quantity / (i^g|V x A^p/2. In particular, we notice 
that the dynamics of the radiated field is 

A = {i7,A} = ^ = Y, 

which is equivalent to the equation for the radiated electric field E^: 

dA 



'^-^^-^ /f + (p - A^ - A)2 + Jd' 



H= c/^gc/WmomVl + (p - A^ - A)2 + / d''q 



E. 



dt 



For a wave co-propagating with the electrons in the z-direction, one can define the fc-mode 
of the wave as follows: 

Afc(q^) = ^ I dze-''^ A{ci), 



where q_L = (x, y) and L is the length of the cavity where the interaction takes place. 
This gives the Fourier expansion in the propagation direction Y(q) = Ylk^k{^.±)^^^^ and 
A(q) = Afc(q_|_)e*'^^. Furthermore, since 

^Afc(qj_) 1 -ikzj:(„ / N 

which is obtained from the definition and linearity of the functional derivative, it follows 
from the bracket ([9]) that 

{Yfc(q_L), A„fc/(q'_L)} = -^4fc'5(q± -ql)- 

Since we also have: 

{Yfc(qx), Y,,(qi)} = {Afc(qx), A,,(qY)} = 0, 



the field part of the bracket turns into 

We now consider the paraxial approximation both for the radiated and external fields, i.e. 
we neglect their spatial variations in the x and y directions, so that they are homogeneous 
in the transverse section S of interaction, and null outside it. We notice that this paraxial 
approximation is the strongest approximation involved in the derivation process. The di- 
mensional reduction crucially depends on it. In addition to this approximation, we restrict 
the derivation to a monochromatic wave, i.e. we only take into account one Fourier mode k 
in the propagation direction, and consider the case of a circularly polarized radiated wave. 
Other modes can be included in the derivation in a very similar way, but we have only kept 
one mode for the sake of clarity of the derivation. These two approximations allow us to 
define the complex amplitude of the wave a such that 

A = — ^ - a*e-'''%*] , 

Y = — ^ [ae'^'e + a*e-'^'e*] , (10) 
v2 

with e = (x + iy) / Conversely, a can be defined as 



where V is the volume of the interaction domain, i.e. V — LS with the above notations. 
Then, since 



Sa 



SYk{q±) kS 
Sa i 



we get 



5A,(qx) S 



so that a and a* are the new conjugate variables describing the radiated field. Hence, the 
bracket turns into: 



^ // 



d SF d SG d SF d SG 



d qd pfmom c r '-irj' 

i fdFdG dF dG^ 



'^'kV V (5a da* da* da ^ 

As for the Hamiltonian, since we have V x A = kA and |Y|^ = A;^|A|^ = k^aa* for the 
vector potential A, the energy of the radiated wave now reads: 



/ 



|YP + IV X A| 



= k^Vaa\ 



where we have used the relations e • e = and e • e* = 1. So that the Hamiltonian becomes: 
jj d^qd^pU 

+|A^|2-2pi-(A^ + A) 



1 + p^ + aa* - iV2{ae''''e - a^e-'^^'e*) ■ A, 

1/2 



+k^Va*a - 



ikS 
72 



dz{ae'''% - a* 6-'''%*) ■ (V x A^). 



B. Reduction to a one-dimensional model 

Here we assume that the external field A^ created by the undulator only depends on the 
longitudinal variable z. If the beam of electrons has been injected in a proper way, we show 
below that the motion is exactly described by a one-dimensional Hamiltonian, allowing for 
a reduced - but exact - description of the dynamics. 
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1 + +aa* - i\f2[a(^^% - a*e-*^^e*) ■ 



This reduction follows from the properties of the Liouville operator 7^ = {if, ■}. Recalling 
that: 

m _ 

^ fmom 

the Liouville operator reads 



|A^|2-2px-(A^ + A) 



1/2 



fn 



+Pz 



d 6 

dz6fnic 



+ 



kV 



OH d 



(A + A^)2 
- - f A + A, 
dH d 



{P± - A, 

. d 



A) 



d 



dz 



A + A„ 



5_ 

dpz 5/moii 



^ da da* da* da 

The Q^-fj^ — term has disappeared since its factor d{A + A^)/9q_L vanishes, as the fields 
are assumed not to depend on the transverse direction q_|_. As a consequence, there is no 
evolution for the particles distribution /mom(q, p) along the direction p^. This can be seen 
by considering a distribution function of the following form 

/mom(q,p) = /(q,Pz)5(p±). 

Under the Liouville operator, we see that 7i/mom is also proportional to S{p±), using an 
integration by parts and the fact that 

<^/mom(q)P) r-/ /w/ /\ 

— —— = 5 q-q <5 p-p . 

'3/mom(q,P') 

We recall that p^ has been translated by A^ in Sec. IIII At so in other words, if the beam is 
initially injected with the transverse velocity A^ 



'1 + |A^p [from Eq. ([Tj)], it remains with 
this specific transverse velocity. Then, once restricted to the S{p±) distribution, it comes 
that the transverse profile of / does not act on the longitudinal dynamics any more. In 
other words, the set JF of observables F which do not depend on the transverse component 
of the distribution, i.e. F such that 

d 6F 

0, 



is stable by Ti. (i.e. TiF G if F G JF), since A and A^ do not depend on q^. This 
allows one to focus on the longitudinal dynamics by defining a reduced dynamics on JF, with 



f{z,pz) the longitudinal distribution as new variable associated with the particles. The 
bracket reduces to : 



dzdpj 



d 5F d 5G d 5F d 5G' 



+- 



I / dFdG dF dC 



kV \ da da* da* da 



and the Hamiltonian becomes 



H 



dzdpzfy 1 + p1 + aa* — iv^(ae*^^e — a*e *^^e*) ■ + | Att,p 
+k^Vaa* J dz{ae'''% - a*e-''''e*) ■ (V x A^). 



Since the motion is now one-dimensional, we drop the label z of the momentum in what 
follows. 



C. Particles-field phase frame 

While Hamiltonian ffTTl) is reduced to distribution functions with one dimension (one in 
space and one in momentum), it still contains some terms which are not specific to the 
interaction like the term aa*. The emphasis can be put on the interaction between the 
particles and the wave by considering the dynamics into the particles-field phase frame. We 
consider a specific medium for the interaction between the particles and the radiated field. 
For example, we can consider a linear undulator, such as in a Free Electron Laser. The 
transverse field produced by such an undulator reads : 

A^ = {e-'''-% + e'^»"e*) . (12) 

First we neglect the effects of finite size : The last term in Eq. fill I) vanishes since it is a 
sum of terms proportional to / dze'^^^'^^'''^^^ terms. The Hamiltonian becomes 

H = J J dzdpf \Jl +p^ + aa* — m^(ae*('^+'^«')^ — a*e~*(('^+^^)^)) + + k'^Vaa*. 

In this Hamiltonian, the last term simply yields the propagation of the electromagnetic 
wave. Indeed, since {k'^Vaa* ,a} = —ika, it generates an e~*^* factor for a{t), and so it is a 
pure propagation term. This remark calls for a time- dependent change of coordinates. This 
procedure is fairly standard, e.g., in quantum mechanics, and corresponds to the interaction 
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representation 15|. Since time is not a variable for this model, we first need to extend 
phase space to add a new pair of canonically conjugate variables with one being similar to 
time. More precisely, we define the pair of conjugate variables {t,E), such that the new 
Hamiltonian and bracket read 



H,^t[f,a,a*,E,T]=H[f,a,a*] + E, 



{F,G} 



dzdpf 



+ - 



d_5£d_5G _ d_5£d_SG' 
dpjfdzjf ~ dzjfdpjf_ 
OFdG dFdG\ OF dG OF dG 
kV \ da da* da* da I dE dr dr dE' 



so that f = {iJcxtjT} = 1, which means that, practically, r is identical to the evolution 
variable t. We now consider the canonical change of variables (/, a, a*, r, E) (/, a, a*, f, E) 
such that 



a = ae 



ikr 



E = E + k^Vaa* 



T = T. 



It can be checked that the symplectic form ikVda* A da + dE A dr is conserved by this 
transformation. This results in a conservation of the Poisson bracket, while the Hamiltonian 
now reads: 



H 



Finally, the dynamics can be studied in the particles-field phase frame, the latter phase 
being defined as ^ = {k + ku,)z — kr, by considering the change of variables 

mp) = f{z,p)/{k+k^), 

a = a. 



E = E + 



T = 



dOdpfp, 
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where / has been divided hj k + ky^ for normahzation purposes. We notice that the term 
—k/ {k + ky,) J j dOdpfp has been added to E in order to ensure the canonicity of the change 
of coordinates which translates here into the condition {E, /} = 0. Here we have used the 
following properties of the functional derivative: If we perform a change of variables of /, 
denoted f = f o g, then the functional derivative of F[f] = F[f] satisfies 

6F 5F , , 

which is obtained in a straightforward way from the definition of the functional derivative. 

Now that we have performed the time- dependent change of coordinates and that the 
Hamiltonian is time- independent, the (r, E) variables are somehow artificial and decoupled 
from the other ones. It is more convenient to work in the reduced space and drop this 
additional pair of variables. For notation purposes, we drop the bars over the other variables. 
The Hamiltonian now reads: 

k 



H = JJ dedpf 
The resulting Poisson bracket writes 

{F,G} = {k + ky,)jj d9dpf 



a/1 + — ia^(ae*^ — a*e + aa* 



-p 



(13) 



+- 



d_S£d_SG _ d_5£d_SG 
d^JfdeJJ ~ mlfd^Jf 
dF dG dF dG 



kV \ da da* da* da 



(14) 



D. Resonance condition and high-gain amplification 

The next step is to expand Hamiltonian fll3p around a resonant value. Following Eq. 
Hamiltonian ( |T3l) and bracket ( |T^ yield the following equations of motion: 




p 



+ — iau,[ae:'-^ — a*e~''-*^) + aa* 



-ie\ 




ay,{k + ky,){ae'^ + a*e-'^) 



df 



da 
It 



1 + p + — iaw{ae^^ — a*e 

'v 



+ aa* dp 

ie 



dOdpf 



'1 + p^ + alj — iaw{ae'^^ — a*e *^) + aa* 

From these equations it can be seen that the system is at equilibrium for a = and f{6,p) = 
6{p — pfi)F{6), where F{6) is a distribution which satisfies J d6e~'^^F{6) = 0, and pr is given 
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by 



PR k Q_ 



This resonant momentum can be linked to a resonant energy 7/j for the particles, defined 
by 



7.^Vl + a^+pi^v'l + <^J=p=. (15) 



In the limit k ^ kw, the definition (jlSll of the resonant energy yields the usual definition 
'jR = y/k{l + a^)/(2/c^) (see for example I2]). 

However, this equilibrium is unstable, and exposed to small perturbations, the wave starts 
growing and destabilizes the particles aX p = pr. This instability is responsible for the high- 
gain growth of the wave, which is taken advantage of in devices such as FEL. The dynamics 
can be linearized around this equilibrium point: Assuming the momenta of the particles 
remain close from the resonant one pr, we shift p by pr^ by defining f{9,p) = f{9,p) with 
p = p — pr. We also consider that the amplitude of the radiated field is weak compared to 
the resonant energy 

\a\ < IR- 

Then, Hamiltonian 0131) expands, at the first order in a and second order in p (we have 
dropped the hat over p and /) , as follows: 



i/u„ = / / dOdpf 



1 + ia. 

apj 

associated with the bracket (fT^ . 



Yr 2 2^R 



The equations of motion can be normalized through the following change of variables 

no\p') = -^f{e = e\p = p'/p), 

a' = ea. 

Moreover, we include a rescaling of time t' = at, i.e. H' = H/a, and we consider a new 
Hamiltonian uH' with a new Poisson bracket z/~^{., .} (which does not change the dynamics). 
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Using 



a 



IR 



2W 



1/3 



^_ 2 f 2kVkl{k + K/2f 
k -\- k^ \ 



2\ 1/3 



4eV^k^{k + k^/2) 



1/3 



the Hamiltonian and the bracket become : 

^2 



H 



jj dedpf--i J J dedpf 



(16) 



dOdpf 



d_S£d_5G _ d_S£d_5G 
d^JJdeJf ~ deJfd^Jf 

Finally, the canonical change of variables (a, a*) — > (0. 



dF dG dF dG 



da da* da* da 
I) such that a = vTe"""^ (so that 
{4>,I} = 1) allows one to retrieve Hamiltonian ([S]) associated with the bracket (jl]). 



E. Conserved quantities 

Apart from H as given by Eq. ([5]), the Hamiltonian system of charged particles interacting 
self-consistently with electromagnetic fields has its total momentum as conserved quantity, 



as it was reported in Ref. 



P[f, E, B] = jj d^qd^pfp + j d^qE x B. 



We perform the same approximations and reductions as the ones done on Hamiltonian (I5l). 
In this way, we recover the conserved quantity of the reduced model. For example, in the 
canonical formulation (IHI) and (IHI), the total momentum turns into: 



P[/^om, A, Y] = J J d'qd'pUoUp - A) - y d'qY x (V x A) . 

Furthermore, when considering a monochromatic, circularly polarised plane-wave such 
as the one given by Eq. fllOl) . the total momentum is decomposed into a transverse and a 
longitudinal component: 

Px[/, jj d'qd'pf ipi_ - A - A^), 

P,[/, A, Y] = jj d^qd'pfp, + eVaa* + k j d'q (A^ x Y) ■ e,. 
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When focusing on the longitudinal dynamics (see Sec. IIII B\} . the conserved quantities 
can be dropped since they provide information on the transverse dynamics. Then, when 
considering the specific external field (fT2l) . if we neglect the finite-size effect, the last term 
of the longitudinal momentum can be dropped, so that we get, for Hamiltonian ( |T3l) . the 
following conserved quantity 

P,[f, a, a*]= jj d^qd^pfp + k'^Vaa*. 

Using the same normalization as in Sec. IIIIDl becomes a conserved quantity for Hamil- 
tonian (fTBjl : 

P = Jj dOdpfp + aa*. 
Finally, using the canonical change of variables (a, a*) y-^ {(p,!), it becomes: 



P = JJ dOdpfp + I, 
i.e. the average momentum plus the intensity is conserved by Hamiltonian ([3]). 

Conclusion 

In this paper, we derived a reduced Hamiltonian for the interaction between a wave and a 
beam of charged particles driven by an external field, on the sole assumptions of transverse 
fields and on-axis injection for the particles. A resonance condition - around an unstable 
equilibrium point - has been identified, allowing for a linearization of the dynamics. Finally, 
under the extra hypotheses of high-energy particles and weak radiated field, a paradigmatic 
Hamiltonian has been retrieved within a fully Hamiltonian treatment. The main advantage 
of the present derivation is a fully algebraic framework which is well suited to include ad- 
ditional effects, like for instance, higher order terms in the expansions, or strategies using 



other harmonics of the radiated field 16|. We have shown here that this treatment allows 
one to recover some general features of the Vlasov-Maxwell equations in a very natural way, 
like a conserved quantity of the flow. 
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